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STABLE ESTIMATES FOR SOURCE SOLUTION OF CRITICAL 
FRACTAL BURGERS EQUATION 


TOMASZ JAKUBOWSKI AND GRZEGORZ SERAFIN 


Abstract. In this paper, we provide two-sided estimates for the source solution of d- 
dimensional critical fractal Burgers equation wt — -|-6• V ('u|m|'^) = 0, q = {a— l)/d, 

a G (1, 2), b G by the density function of the isotropic a-stable process. 


1. Introduction 


Let d eN and a G (1,2). We consider the following pseudo-differential equation 

f Ut — + b ■'V {u\u\'^) = 0, t > 0, X G /.. .., 

\ u{0,x) = M6o{x), 

where M > 0 is arbitrary constant and 6 G is a constant vector. In this paper, we 
focus on the critical case q = {a — l)/d. Here, denotes the fractional Laplacian 

dehned by the Fourier transform 

= -Krm, <!> 6 cr(iRb- 

Equation fll.ip for various values of q and initial conditions uq was recently intensely 
studied ([asiEiiz]). For d = 1, the case g = 2 is of particular interest (see e.g. 

[131 [H]) because it is a natural counterpart of the classical Burgers equation. Another 
interesting value of q is In [1] authors proved that the solution of fll.llh which we 
denote throughout the paper by UM(t,x), exists and is unique and positive. It belongs 
also to (M'^) for every p E [1, cxo]. The exponent q = is critical in some sense. It 
is the only value for which the function UM{t,x) is self-similar. Is satishes the following 
scaling condition ([!]) 

UMit, x) = ax), for all a > 0. (1.2) 

Furthermore, the linear and the nonlinear terms in fll.ip have equivalent influence on the 
asymptotic behavior of the solution. If g > (a — l)/d, the operator plays the main 
role. More precisely, for such g and a function u satisfying fll.ll) . with not necessarily the 
same initial condition, we have 


lim 

t—)-oo 


u{t, ■) — ■) 


0 , 


for each p E [1, cxo]. 


For g < {a — l)/d another asymptotic behavior is expected. In addition, taking g = 
for d = 1 and a = 2 we obtain the classical case, which makes the equation fll.ll) with 
critical exponent g one of the natural generalizations of the Burgers equation. 
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Till the end of the paper we assume that d > 1, a G (1, 2) and q = Let p(t, x) be 
the fundamental solution of 

Vt = (1.3) 

In [7] the authors proved that for sufficiently small M there is a constant C = C{d, a, M, b) 
such that 

UMit,x) < Cp{t,x), t > 0, a; G (1.4) 

In this paper we get rid of the smallness assumption of M. Furthermore, we also obtain 
the lower bounds of um- We propose a new method which allows us to show pointwise 
estimates of solutions to the nonlinear problem fll.ip without the smallness assumption 
imposed on M. This method has been inspired by the proof of [HI Theorem 1]. Our main 
result is 

Theorem 1.1. Let d > 1 and a G (1,2). Let UM(t,x) be the solution of the equation 
M.l\) with q = There exists a constant C = C{d,a, M, b) such that 

C~^p{t,x) < UMit^x) < Cp{t,x), t > 0, X G M'^. 

The fractional Laplacian plays also a very important role in the probability theory as a 
generator of the so called isotropic stable process. The theory of its linear perturbations 
has been recently significantly developed, see e.g., laiHiiiniiiiiiiTiiiiiiisiiHiE]. However, 
since the term b- V{\u\'^u) in fll.ip represents a nonlinear drift, methods used in the linear 
case often cannot be adapted. In the proofs we mostly use the Duhamel formula and its 
suitable iteration. The scaling condition fll.2p is also intensively exploit. 

The paper is organized as follows. In Preliminaries we collect some basic properties of 
the function p{t, x) and introduce the Duhamel formula as well. In Section 3 we prove 
that the solution of fll.ip converges to 0 as |x| —?■ oo. In section 4 we prove Theorem 11.11 


2. Preliminaries 


2.1. Notation. For two positive functions f,g we denote f ^ g whenever there exists a 
constant c > 1 such that /(x) < cg{x) for every argument x. If / ^ and g ^ f write 
f ^ g. If value of a constant in estimates is relevant, we denote it by fc G N, and it 
does not change throughout the paper. 


2.2. Properties of p{t, x). The fundamental solution of fll.3l) may be given by the inverse 
Fourier transform 

p{t, x) = (27r)-'^ f t > 0 , x G . 


This implies the following scaling property 

p(t,x) = a'^p{a°‘t,ax), for all a > 0. (2.1) 

Note that mm(L possesses exactly the same property. Let p(t, x, y) := p(t, y—x). Below, 
we give two well-known estimates of p and the gradient of p (see [H] for more details). 


p{t,x,y) ^ 

Vyp{t,x,y)\ ^ 


t 


(fV“ + \y 

t\y 


-x|)"+“’ 
— x| 


((!/«+ 


( 2 . 2 ) 

(2.3) 


We will need the following lemma 
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Lemma 2.1. Fort,e > 0 we have 


'B{0,e) 


Proof. By formula fl2.2p . 


'5(0,e) 


p{t, 0, w)dw 


tdw 


'5(0,e) 

re tr‘^~^dr 


= c 


'o (tv« + r) 


d+a 


t 


'1 ^ d -1 


r“ 


= c ■ 


^“./o /"'Md + r) 


d+a ' 


(2.4) 


It t'/” > £. we estimate the denominator in the integral by and we get 


'5(0,e) 


p(t, 0, ss 


(2.5) 


In the case t^^°‘ < e we substitute r = in fl2.4p . then 


'5(0,e) 


p(t, 0, 


~ / - ~ 


'0 (1 + m ) 


d+a 


( 2 , 6 ) 


Combining fl2.5p and fl2.6p . we obtain the assertion of the lemma. 


□ 


2.3. Duhamel formula. Our mail tool is the following Duhamel formula, 


UM(t,x) = Mp(t,x) + 



p{t — s, X, z)b ■ Vz[um{s, dz ds. 


(2.7) 


In the following, we assume that UM{t,x) = is a nonnegative self¬ 

similar solution of the equation (12.71) such that um{^, ■) E IP (Mf) for each p G [1, cxo]. As 
it was mentioned in Introduction, the existence of such a function was shown in [1]. 

It turns out that the integral in (12.7p is not absolutely convergent, but integrating by 
parts we obtain a more convenient form 


UMit, x) = Mp(t, x) 



b ■ Vzp{t — s, X, z) [um{s, dz ds, 


( 2 . 8 ) 
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which is absolutely convergent. Indeed, we have \\um{s, ■)\\p < s p g [l,oo] 

(see i)- Hence, by fl2.2p and fl2.3p . 

ft f 

\b ■ V zp{t — s,x, z)[um{s, dzds 


{V zPit — s, X, z)\[uMis, z)]'^^^ dz ds + / / zPit — s, x, z)\[uMis, z)Y~^^ dz ds 



>0 Js.‘i 
rtl2 


< 

r\j 


i/2 -'IR'* 

{t - s)-^/^p{t - S, X, 2)s-'='(9+b/« dz ds 


/•i/2 /• pt 

i / {t- s)-^^+^'>/‘^[uM{s,z)Y+^dzds + 

0 jR‘t Jt/2 jRd 

pt/2 pt 

</-(d+l)/a / ^-(a-l)/a^^ ^/-%+!)/« / -1/a < /-d/a_ 

io Jt/2 

Now, using the scaling property fll.21) . we obtain 

b ■ VzP{t — s, X, z)[um{s, zyf^'^ dzds 

VzP{t - s, X, dz ds 

^(l-Q:)/Qy^^j'/. _ dw ds 




'0 



'0 JRd, 


= a 


'^wP(t — r", X, rw) dw dr. 


'0 


Finally, we get 


r-il/° 


tiM(^, x) = Mp{t, x) — a 


'0 Jr‘^ 


— r“,x,rt(;)[tiM(l,'«^)]^''’^c^'?Fdr (2.9) 


3 . Behaviour of um{ 1 ,x) at infinity 
Due to the scaling property fll.2p it suffices to consider UM(t,x) only for t = 1. 
Lemma 3.1. There is a constant Ci > 0 such that for every x G 'Mf, 

Um{^,x) < Mp{l,x) + Cl [ [ {1 — r°')~^^°'p{l — r°',x,rw)[uMi^,w)]‘^~^^ dwdr. (3.1) 


'0 


Proof Formulae fl2.2p and 02.31) imply that 

|V2p(l — r, X, 2;)| < (1 — r)“^'^"p(l — r, X, 2 ;), rG [0,1), x, 2 ; G M'^, 
and the assertion follows by 02.91) . 

Now, we will show that the function um{^,x) vanishes at infinity. 
Proposition 3.2. IFe have lim um{)-iX) = 0. 

|ir|—>-oo 

Proof. Let us rewrite 03.ip into the form 

“m(1, x) < p(l,x) + /i(x) + hix), 


□ 
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where 


Ii = 


h = 


.1/2 


' 1 / 2 , 


(1 — r“) ,x,rw)[uM{^,w)Y^^ dw dr, 

(1 — — r“, X, rw) dw dr. 


Note that qr := w)Y^'^dw —)■ 0 as i? —)■ cxo. Using fl2.2p and estimating 

1 — r“ PS 1, r G (0,1/2), we obtain for |x| > R, 


»l/2 


Il = 


.. .dw dr + 


< 

r\j 


'0 JB(o,n) 
/.1/2 /. 


'o Jb{o,r) 


»l/2 


'o Jb{o,r)'^ 


... dw dr 


X 


d-\-a 


mm(1, dw dr + 


.1/2 


'o Jb(q,ry 


[^^(l, w)Y^^dw dr 




\RY+<^ 

which is arbitrary small for sufficiently large R. Consider now the integral R and fix 
e > 0. There exist -Ri, -R 2 >0 such that 


lB{0,RiY 


p{l, 0, w)dw < e, 


' B{0^R2)^ 


mm(1, w)dw < e 


d+l 


The latter inequality implies that the measure of the set {w G 5 ( 0 ,-R2)'^ : um{^,w) > e} 
is less or equal to It gives us for r G (1/2,1) and |x| > Ri + R 2 

/ p{l — r^,x,rw)[uM{^iW)Y^'^ dw 


p{l — r^", x,w)[um{^,w/ dw 


< 


|■'IM(1, 0 


|l+<? 

I CXD 


p(l — r", X, w) dw + 


-l+g 


I|wm(i,-)I1 


' B{x,RiY 

1+9 f 

°° J B{x,Ri), 

UMY,w/r)>6 


I B{x,Ri), 
UM(}-,w/x)<e 


p{l — r“, X, w) dw 


p{l — r“, X, w) dw. 


Using scaling property fl2.ip and substituting tc = x + (1 — we get 




p(l — r“, X, w) dw = 


's(o,iii(i-r“)-i/ay 


p(l, 0, z) dz < e. 


|{m; G B{x,Ri) : UM{l,w/r) > e}| = 

= |{u; G B{x/r, Ri/r) : Um0-,w) > e}| 

< |{m; G 5(0,52)'' : mm( 1 ,mi) > e}| 

<e'^ < |5(0,e)|. 


Moreover, 
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£ 

(1 - + £ 
-l/a I ^ 


Hence, by the fact that p{s,x,w) is a decreasing function of \x — w\ and by Lemma I^TTl 
we obtain 

/ 2 (x) < ^ + / / p{l — r'^,x,w) dwdr 

J1/2 JB(x,e) 

<e+ ^ dr 

dl/2 

rl-{l/2r 

<e+ - ) du 

Jo \uy» + ej 

= e+ / ...du+ / ...du 

Jo Je “/2 

pd/2 /•l-(l/2)“ 

<e+ u~^^°‘du + 7 -: / u~^^^du 

~ Jo (1 + 

<P^pi»-m^^d/2_ 

Therefore, for sufficiently large \x\, integrals Ji and I 2 are arbitrary small. This ends the 
proof. □ 

4. Proof of Theorem 11.11 

In this section, we prove the main theorem of the paper. First, we dehne some auxiliary 
functions. Let 

H{x,w)= f (1 — r")“^/“p(l — r", X, rtc)(ir. (4.1) 

Jo 

For (3 G (0,1), we dehne 


H{x,w) = / r ^(1—r") ^^°‘p(l — r°‘,x,rw)dr. 

Jo 


Additionally, for i? > 0, we denote 

huix) = / H{x,w)[uM{l,w)/^'^dw, 

Jb(o,r) 

hR{x)= H{x,w)[uM{l,w)/^‘^dw, 

Jb{o,r) 

Hr{x)= / H{x,w)uM{^,w)dw. 
Jb(o,ry 

Note that H{x,w) < H{x,w) and hji{x) < hji{x). 

Lemma 4.1. For x G we have 

/ H{x,w)p{l,w)dw = C2 p{3,x), 


(4.2) 

(4.3) 

(4.4) 

(4.5) 


(4.6) 




TT 


a sin (tt/q!) 


> 1. 


where 
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Proof. By scaling property and Chapman-Kolmogorov equation for the function p{t, x, y), 
we get 

/ p{l — s°‘,x, sw)p{l,w) dz = / s~'^p{s~‘^ — 1, s~^x,w)p{l,w) dz 

= s-‘^p{s-‘^,s-^x,) = p{l,x). (4.7) 

Consequently, 

/ H{x, w)p{l, w)dw = 



(1 — s") ^^°‘p{l — s°',x,sw)p{l,w)dwds 

= p(i, x) (1 - s")"^/"ds = ir r p(l, x) 

TT 

p(l,x), 


a sin (tt/q;) 

where the last equality results from the Euler’s reflection formula. □ 

The next step is to provide a Chapman-Kolmogorov-like inequality involving functions 
H{x,w) and H{x,w). At hrst, we present a technical lemma. 

Lemma 4.2. Let ft > 0 be fixed. For v E (0,1), we have 

f _ y^yP»dr ^ v-^{l - 

J V 

Proof. Denote the above integral by I{v). Since — W ^ {a — b)a'^~^ for a > 6 > 0 and 
7 > 0 (cf. Lemma 4 in [IS]), we have 1 — r" ps 1 — r and r“ — (r — T)r““^. Hence, 

I{v) ^ f - r)-^/“(r - v)-^^^dr. 


For V > 1/4, we estimate ^ ^ 1 and substitute r = 1 — u{l — v), which gives 


us 


I{v) ~ (1 — t) 




u 




In the case v < 1/4, we split the integral into and obtain 

/• 1/2 


I{v) 


= V 


ri/“-i-/^(r - T)-i/“dr + [ (1 - r)-^/"dr 

Jl/2 


o 


a — 1 


V ^ + 1 T 


which is equivalent to the required formula. 

Since a > 1, we immediately obtain the following 
Corollary 4.3. Let /3 > 0 6e fixed. For v E (0,1), we have 

[ r-^(l - r“)-^/“(r" - i;")-^/“dr < ^-^(l - w)-^/“. 


□ 


This allows us to prove the following lemma. 
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Lemma 4.4. There exists a constant Cs > 0 such that for x,z & we have 

f H{x,w)H{w,z)dw < C^H{x,z). (4.8) 

Proof. Scaling property and Chapman-Kolmogorov equation for the function p{t, x, y) 
give us 

[ H {x, w)H{w, z)dw 

1 — sw){l — — r°‘,w,rz)dw ds dr 

— 1, — r^,w,rz)dw ds dr 

ayl/as-dpf^S-'^ _ ^ 

r-^(l - - s"r", x, srz)ds dr. 

Substituting s = v/r in the inner integral and then using Fubini-Tonelli theorem, we get 

/ ^ {x, w)H{w, z)dw 
JM.d- 

= [ p{l-v‘^,x,vz) [ r-^(l 



1 pi 



0 Jo 


By Corollary 14.31 we obtain 04.81) . 


□ 


Remark 1. Lemmawhich plays an important role in the above-given proof of Lemma 
14.41 does not hold for (3 = 0. This explains partly the form of the functions H{x,w) and 
hnix). 

As a consequence, we get 

Corollary 4.5. For x G we have 

[ H{x,w)hR{w)dw < Cshnix), (4.9) 

jRd. 

[ H{x,w)HR{w)dw < CsHnix). (4.10) 

Now, we pass to the proof of the main result of the paper. 

Proof of Theorem II.11 Let Co = Ci(C 2 VC 3 ). By Lemmawe may choose p G (0, 1) 
and i? > 0 such that \um{'^,x)\ < ( ^ ) 1^1 > Thus, by Lemma IXTl we have 

CiV f 


mm(1 , x) < Mp(l, x) CihR{x) + 


C 


0 Jb(o,ry 

< Mp(l, x) + CJir{x) + ^ Hr{x). 

Go 


H{x,w)uM{^,'w)dw (4.11) 


(4.12) 
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We put fl4.12p to fl4.1ip and, by Lemma 14.11 and Corollary 14.51 we have 

um{1,x) < Mp{l,x) + CihR{x) (4.13) 


CiV 


0 Jb(o,ry 


H{x, w) 


Mp{l, w) + CihR{w) + ^^Hr{w) 


dw 


< M(1 +T])p{l,x) + Ci{l + p)hR{x) +p‘^Hr{x) . (4.14) 

Now, we put fl4.14p into fl4.1ip and, by Lemma 14.11 and Corrolary 14.51 we obtain 
mm(1, x) < M{1 + p + p‘^)p{l, x) + Ci{l + p + p^)hR{x) + p^Hr{x) . 

Hence, by induction. 


mm(1, x) < h^p(l, x) + Ci'^ r]^hR{x) + p'^^^HRix ). 

Taking n —)■ cx), we get 


fc =0 


k=0 


, , M . , Cl ~ . 

um{1,x) < - - p{l,x) + -- hR{x). 

1 — p 1 — p 


(4.15) 


Furthermore, since both functions p{l, •) and mm( 1, •) are continuous and nonnegative 
(see [1], proof of Theorem 2.1), they are comparable on every compact set. Hence, we 
focus only on large values of |a;|. 

We first prove the upper estimate. Let |a;| > 2R. For \w\ < R and s G (0,1), we have 


\x — stcl > |x|/ 2 , and consequently p{s,x, sw) < |x 


|™|—d—o 


p{l,x). Hence, 


n 


fiRix) = 


Jb{o,r) Jo 

<p(l,x) ||nM(l,-)ll^'^ 


r ^(1—r“) ^^°‘p{l — r°‘,x,rw)[uM{^,uj)Y~^'^dr dw 
f r~^{l — r°‘)~^^°‘dr dw 


Jb{o,r) Jo 

~ ll'«M(l,-)lli^‘' P{^,x), 
which, by fl4.15p . gives us 

um{1,x) <p{l,x). (4.16) 

To provide lower estimates we show that the integral in 02.911 is much smaller than Mp{t, x) 
for sufficiently large |x|. By 04.16p and 02.311 . there are constants C 4 , C 5 > 0 such that 


a b ■ V.ujp(1 — s", X, sw) [mm(1, w)]'^’'"^ dw ds 



< C 4 


<C5 



I V^p(l — s", X, stc)[p(l, ta)]^’'“^| dw ds 
1 

lo JB{o,R) (1 - + \x-sw 



p{l — s", X, sw)p{l, w) dw ds 


+ Cs 



(1 — s“) ^/“p(l — s“, X, sw) [p{l, dw ds 


0 J B(0,RY 
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for any ^ > 0. We take R such that p(l, w) < (1 — s“) for \w\ > R. 

Then, for |a:| > {2R) V ( 8 C 5 /M), we obtain 


a / b ■ \/wp{l — 3°", X, sw)[um{^,w)]'^~^^ dw ds 


0 JM.d, 


< 


2C, 





p{l — s“, x, sw)p{l, w) dw ds 


'0 Jb{o,r) 

M 


4 - s“)-V«rfs Jo 

Hence, by fl4.7p . we get 




(1 - s“) 


'R(0,R)= 


p(l — s°‘,x, sw)p{l, w) dw ds. 


a / b ■ VviPi^ — s°‘,x, sw)[um{^,w)Y~^^ dw ds 


IQ 

which, together with fl2.9p . implies 


M , , 


M 


uni^^x) > —p{t,x). 


The proof is complete. □ 

Remark 2. Using the result of Theorem M. 1\ as well as the method of proving it one can 
show that there exists a constant C = C{d, a, M, b) such that 

\VxUMit,x)\ < Ct~^/^p{t,x). 


However, this estimate does not seem to he optimal. 
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